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Abstract 


An extension of the method of moments is developed for the numerical integration of the kinetic equa- 
tions of droplet spectra evolution by condensation/evaporation ^jid by coalescence/breakup processes. The nu- 
mber density function ni.(x,t) in each separate droplet packet between droplet mass grid points .) is 

represented an expansion in orthogonal polynomials with a given weighting function. In this way droplet num- 
ber concentrations, liquid water contents and other moments in each droplet packet are conserved and the prob- 
lem of solving the kinetic equations is replaced by one of solving a set of coupled differential equations 
for the number density function moments. The method is tested against existing analytic solutions of the cor- 
responding kinetic equations. Numerical results are obtained for different coalescence/breakup and condensati- 
on/evaporation kernels and for different initial droplet spectra. Also droplet mass grid Intervals, weighting 
functions and time steps are varied. 


Introduction 


There are three major difficulties in numerical computation of droplet spectra evolution by condensa- 
tion/evaporation and coalescence/breakup processes, which occur as a result of interaction between droplets 
and vapour-air environment. 

a) The relaxation time of condensation/evaporation process is much smaller than the relaxation time of 
coalescence/breakup process and therefore the numerical computations of these processes require a different 
time steps, 

b) A correct approximation of droplet number density function nj^(x,t) is required in each separate 

droplet packet between d’-oplet mass grid points in Berr^ s approximation^*^ of number density 

function neither the number concentration nor the liquid water content of droplets are conserved. In Block's 
method^*^ it is impossible to estimate the error of the approximate numerical solution and an assumption 
is made *i .t in sach separate droplet packet all droplets are spread over the whole mass interval 
This Bleck's uniform distribution hypothesis gives as a result a significant increase of the mass conver- 
sion velocity from small droplets to large drops. 

c) The problem of correct cocputation of water vapour supersaturation taking into account the release 
of latent heat of condensation/evaporation during the time step used for the numerical computation of the 
condensation/ evaporation processes . 


Method of mosients 

In this study an e.xtension of Block's method^ and of the method of moments ^*^ is developed for nu- 
merical Integrating the kinetic equations of droplet spectra evolution by condensation/ evaporation and 
coalescence/breakup processes. Each separate droplet mass interval .where Xi.^iasxy is considered 

as a droplet packet with its own number concentration, liquid water content^and other mraents. The unknown 
number density -function n^(x,tj in each droplet packet is represented by an expansion in orthogonal pol)mo- 
mials with a given weighting function 


"kCz.tj-W^Cz.tJI^ai^CtjG.^d) (1) 

where z>x/Xv represents the nondimensional mass of droplets in the packet (XviXt..,). 

Wjj(z,t)--weighting function. 

^ik(z)-are pol)momials orthogonal in the range (l,s) with weighting function N|^(z,t). 

aij^j-tj-are the expansion coefficients which are expressed as linear combina..lons of the number 
density function moments and depend on the moatents of Nj^(z,t) as well. 

In this way droplets number concentration, liquid water content and other aximents are conserved in each 
separate droplet packet and the problem of solving the kinetic equations is replaced by one of solving an 
infinitive set of coupled differential equations for the number density function moments. However, approxlsM- 
ting ng(z,t) by means of first L terms of the expansion (1) and also replacing the expansion coefficients 
aik(t) by means of linear comblnttlons of number density function moments we obtain a finite set of 
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coupled differential equations to conpute the first L aioaients of the n, (i,t)^. Note that the approxlaa- 
tlon of the njj(z,t) by neans of the first L terms of the expansion (1)*^ assumes that the expansion 
coefficients for 1 i., iaL«l, are zero, that Is 


( 2 ) 

By (2) the larger order moments of the number density function may be determined by means of the 
first L moments of the nj^Cz.t). 

The accuracy of the approximation of nj^Cz.t) by means of the first L terms of the 

e^^ansion (1) in droplet packet( depends on the choice of W. (z,t) as well as on the 

choice of droplet mass grid points. Tne number density function n(x,t) Is a distribution funct- 
ion which in experiments is determined only for the range (x,x*d)0> If droplet mass increment 

dx is much smaller thin x, the n(x,t) in each range (x,x+^) represents a piecewise constant 
function and an arbitrary moment of n(x,t) in the range (x.x«dx) is expressed by means of the 

zero-order moment dN(x,t)>n(x,t)dx; therefore also the number density function n(x,t) is determined 

by zero-order moment d»(x,t) and does not depend on the choice of the weighting functions in 

the range (x,x*dx). 

For the numerical computations of the kinetic equations it is imposs'ble to cho- 
ose droplet mass grid points for which CX]^,.j-x.) is much smaller than Xj^ ; therefore for such 

grid points the zero-order approximation of the^ expansion (1) (L« 1) will ba incorrect for an 

arbitrary weighting function Wj^(zvt). So, if for numerical solution of the corresponding kinetic 
equations we have choosen the droplet mass grid points Xj.^,»sxv , then we have to define for 
what value oi the grid intervals the first-order approximation of (1) (L“2) will be 

correct that is will be not dependent on^ the choice of the weighting function **v(z,t). It is 
shown that for s equal or smaller than two the first-order approximation of the expansion (1) 
that is the approximation of nj^(z,t) by means of the first two moments is sufficiently correct. 

The relations are obtained to compute the zero-order and the first-order approxima- 
tions of the polynomial expansion (1) for an arbitrary range (l,s) and for an arbitrary weight- 

ing fume ions. 

It should be noted that the first-order approximation of the expansion (1) (L>2) 
describes not only the case where the droplets in the packet (Xk,x.^.) are spread over the 

whole mass interval (*k»*v+i) > ®I®° where the droplets*' in the packet (*k»*ic*i) 

are located only in the ''part of the whole mass interval 2) ^ 


Numerical method for integrating the condensation/evaporation 
kinetic equation 

In this study a separate treatment of microphysics of condensation/evaporation process 
(without advectlon phenomena) is adopted and this process for sufficiently small time steps is 
cor.sidered as a space-homogenous process. Such consideration is based on the assumption that 
the tine integration is broken up into separate treatment of the dynamic tendency and of the 
microphyslcal processes^ which control the vapour supersaturation field. 

Differential equations which describe the microphysics of the droplet spectra evo- 
lution by condensation/evaporation processes can be written as : 


9n(x,t) j . 

3t * 77 ' "U.t; J 

(3) 

dx -K - Bxl/^ ♦ tfclx^/^ 

(4) 

S(t)-Q(t)-Q j(t) 

(S) 

Q{t)*I.(t)»Q(0)H(0) 

(6) 
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P,C dT , dMCt; 
r 1 p ~3f *-e^t 


( 7 ) 


dT T 



( 8 ) 


where S(t)-water vapour supersaturation; M(t)> droplets total liquid water content; Q(t)-water vapo- 
ur density; Q (t)- water vapour saturation density at temperature T(t). The equation (3) represents 
the condensation/ evaporation kinetic equation; (4) -represents the individual droplet dlffuslonal growth 
aquation which includes the terms due to surface tension (B) and solute effects -4 t.c} ; (6) re- 
presents the isass conservation equation; (7)- the first law of thenoodynaalcs ; (8)- the Clapeyron 
-Clausius equation. 


Multiplying (3) by *“dx and Integrating from Xj^Ct) to 
account that n(x,t)dx«nQ(y)dy we have the set of number density 


X]^^l(t), and also taking into 
function moments equatiom- 


m 

dM ^ft) 
dt 




ffl-1 

ft / 3^ 




ViW 

(dx/dt)n(x,t)dx« m / UCy.t)) 

Xk(0) 


(dx(y,t)/dtji^ (y)dy 


(9) 


where 

M ®(t)» ; x“ n(x,t)dx ^ / (x(y,t))® ngCyjdy 
x^(0) 

represents the m-order moment of the number density function in the droplet packet with nonfixed 
grid points ( Xjj(t), Xi.*i(t) )> that is in the droplet packet which for t»0 is contained within 
grid points ( Xj^CO), XjJ^j(O) ) »nd for time t transfers in the droplet packet ( Xj^(t), j 

as a result of the droplets dlffuslonal growth. 

For t>0 we have Initial conditions: 
n(x,0j - n^Cy); x(y.0)-y; xitCOJ-x^; *k*lC°)“*k*l 
From equation (9) we hwve for m*0 and m«l 


dN^(t) 

dt 


■ 0 ; 


dt 


/ (d*(y.t)/dt)n^(y)dy 
*k(0) 


CIOJ 


where Nj^(t)«M^(t) and M.(t)eFC (t) represent the number concentration tmd liquid water conte.it 
respectively in the droplet picket with nonfixed grid points (Xj^(t), )• Using (10) we 

have for the droplets total liquid water content 


«ii.i 

« k-1 


Vi(0) 

/ (dx(y,t)/dt)n (y)dy 

Xk(0) 


(U) 


where J-the total number of droplet packets. 


332 



Picard's method of successive approxlmtlons is used for the integration of differen* 
tlal equation (4) with the initial condition for t«0 xoy. It is assumed that the droplet indi- 
vidual diffuslonal growth rate in the (j'»l)-th approximation is determined by the droplet '^ass 
in j-th approximation. Substituting (4) in (11) and approximating nQ(y) by means of the first 
two moments Nj^(O) and M. (0) we obtain the first relationship between total liquid water con- 
tent of droplets and between water vapour supersaturation: 


where 



ill 


dt 


(j) (3-1) (3-1) 

s ( t )* ( t ) ♦ f ( t ) 


.(3-1 


!t) 


J *krl(®) (i-1) 2/3 

I / ((x^J '^(y.t)) ^V(y)dy) / 

k«l X. (0) ® 


(3-1) 1/5 

(X (y,t)) 


(3-1) 

t (t) 


K Z 
k>l 


1/ 1 

(j-1) 1/5 (l-n 2/3 

/ (dyj((-B(x^-’ \y.t)) ♦e(c)(x^^ ^v.t)) 

Xt(0) 




are the condensation/evaporation integrals. 

The second relationship between M^^^(t) and S^^^(t) can be derived using the equations 
(5)-(8). Expanding Qj(t) in Tailor’s series, we have for sufficiently small time steps, ( that is 
for time steps for which ( Tlt)-T(01)« T(0) is satisfied) 


S(3) (t)«S(0)-(l*FjfM^J^(t) -M(0)) (12) 


where 


F- ( ( (Q, (T (0) ) Lg)/ (r jCpT (0) ) ) ( (Lg/RyTiO) ) -1) 

and Lg rmpresents the latent heat of water vaporization. Such an approach gives the possibility to 
obtain the analytic relationships to compute the unknown values of x(y,t), x.(x.(o),t). M]^(>.), 

M(t), S(t), Q(t), T(t), Q (t) for sufficiently small time steps t, (0<t<i.) which*' g,ay be much 
smaller than the time step used fox- the numerical computation the coalescence/brenkup pro- 
cesses. The redistribution of the liquid water mass from the dro,/let packets between nonfij:ed 

.grid points (*) 5 (t) ,X)^+j (t)) among the droplet packets located between fixed grid points (Xjj.Xj^^j) 

which are used for the numerical computation of the coalescence/breakup processes can be computed 
without too much difficulty^. 

Note that according to (10) the number concentration h|^(t) in the droplet packets with 
nonfixed grid points (Xk(t^ iXjj^j (t)) is constant and for initial “^monodisperse droplet spectrum 

M(t) ■ '!i.«))x(t) where N(0) represents the total number concentration of droplets. Therefore for the 
initial monodisperse droplet spectrum and for the case B«0 and 4(c)>0 in equation (4) there 

exists the analytic solution of the set of differential equations (4), (11), (12) which may be 

used for the test of the numerical method developed in this study. 


Numerical method for integrating the coalescence/breakup 
kinetic equation 


The 

can be written 


coalescence/breakup kinetic equation for the droplet 
as: 


number density function n(x,t) 


In 


3n(x,t) 

it 

this 


- --n(x,t) 
equation 


“ x/2 

/o(x.y)n(y,t)dy* /o(x-y,y)n(x-y,t)n(y,t)dy-n(x,t)P. x)» 
*1 *1 

o(x,y) lepresents the coalescence kernel for two 


/n(y,t)Q(y,x)P(y)dy 

X 


(15) 


droplets of mass x and y; 
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i ls th« sMllut BUI in the tpectnia of oropleti', P(x)- the probeblllty thet a droplet of 
St X will break tq> during a unit tiaa and Q(y,x)-tha droplet maaber density function, 'oraed' 
due to breakup of a parent drop of aast y. 

Multiplying (13) by x*dx and integrating (13) froa Xj^ to yields 



\ - »*k k 2 

! x*n(x,t)dx I o(x,y)n(y,t)dy ♦ /x*dx I o(x-y,y)n(x-y,t)n(y,t)dy 
*k ’‘l *k *1 


**k s*^ . 

- / x“n(x,t)P(x)dx ♦ I x*dx / n(y,t)Q(y,x)P(y)dy (14) 

*k *k * 


where 


I. * 

■ * B n B 

M (t)« / X n (x,t)dx ■ X. / 2 n. ( 2 ,t)d 2 
■' *k 1 


(IS) 


represents the B-th order Boaent of the nuaber density function of droplets in the packet 
(X|^,X|^^.). According to Block the first two double integrals on the right tide of (14) can 
be cOBputed sectional ly over snail areas. Using this idea also for the computation of the last 
two integrals in (14) , we have 


at 


J a k*l B k-1 B J B B 

'ii/i.kiil^i.k-l.kii/i.k.k i5k“i.k ‘*k 


(16) 


where k«l, 2, ,J; J is the nuaber of droplet packets; n>0, 1, 2, and b is the order 

of Boaents. In the set of differential e<}uatirns (16) Af^k > *i.k.l,k • ^ k k coalescence 

double integrals which describe the kinetics of the coalescence piocess. Thi 'last two tern* in 

(16) describe the droplet breakup process. 

The computation of the coalescence and breakup double integrals requires the' values 
of the unknown nusd>er density function n(x,t) in each separate droplet packet. Therefore approxi- 
nating again the nuBber density function n,^(x,t) by neens oi the first L teras of the expansi* 
on ( 1 ) and replacing the expansior coeffi tents a>. (t) by aeans of linear conbinations of M|^”(t), 

we have froB (lo) a finite set nf co>ipl<.J differential equations to coeipute the first L ^Seo> 
Bents of the it|^( 2 ,t) in each separate droplet p;.cket. 

There exist conservation relations between the coalescence double integrals, which for 
siaO and si^l describe the droplets nuaber conceritratir.i conservation and liqtiid water lontent con- 
servation respectively during the coalescence interection between two droplet packets. 

The detail;^ Bethod for mnerical integrating of the coalescencw/breaktqp kinetic equ- 
ation is contained in 


Niaerlcal corputetions and co.icluding 
reaarks 


The unknown nuBber density function 0 ^( 2 , t) In each separate droplet packet (^,^ 41 ) 

(Xk 4 i*sXk) is represented by an expansion in orthogonal polynonri’ls with a given weighting func- 
tion in the range (^k'^ktl)* pxobleB of solving the condensation/evaporation and 

coalescence/breakup kinetic equations is replaced by one of solving a set of coupled differential 
aquations for the anaents of the nuBber density function np(s,t). 

The Bethod of Boaents developed in this study, is tested against existing analytic 
solutions of the correspondii ^ kinetic aquations. Nwericsl results *;* obtained for different 
coalescence/breakup kernels, foi 'ilfferert individual droplet diffu*lonal growth rate and for dlf 
ferent Initial droplets spectra. Also droplet Bass gr'.^ points intervals, weighting functions and 
integration tisM steps are varied. 

The results of the nuBerical coBputations of the droplet spectra evolution by condensa- 
tion/evaporation processes Indicate that the convergence of the kicard** Bethod of the succasive 
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q>proxiMtions for the equation (4) for SMll time intervals is sufficiently rapid. For monodis- 
perse initial droplet spectra a comparison betueen nmerical results coaputed by the method of 
this study and between existing analytic solution of droplet spectra evolution by condensation pro- 
cess is made. It is shown that a simultaneous use of the Ficard«s method and of the method of 

moments gives reasonable results which are very close to the existing analytic solution. 

For the mawrical integration of the coipled differential equations for N|^(tl and 
Ng(t) which are derived from C16)> the forward time differencing is adopted. For practical com- 
putations the droplet mass grid points are the most convenient. The results of the nu- 
merical computations indicate that if the chosm weighting function is close to the unknown num- 

ber density function nj^(z,t), the difference between numerical results computed by means of the 
zero-order aporoximation and by means of the first-order approximation of the ejqiansion (1) is 
sufficiently small. The results of numerical co^nitations also show that the difference between 
numerical results computed by means of the first-order approximation of the ejqwnsion (1) with 
different l^(z.t) i^ not large and even if the chosen weighting functions are not close to the 

unknown number d- .-'.ty function n|^(z,t]. the use of the first-order approximation of the expan- 
sion .'1 for di.f'rent Rb(z.t) gives nuawrical results tdtich are close to the existing analytic 

solutions. Therefore one important aspect of the method sf moments is its ability to estimate 

for a ^.<en kernels the efficiency of the choren weighting function, that is the sufficioicy of 
the z. '5-order spproximation of the expansion (1) tdten there is no analytic solution. A prelimi- 
nary numerical results axe obtained to study the influence of the condensation/evqwration proces- 

ses on the ditplet spectra evo’ ition by the coalescence/breakup processes. 

Ibus, the convergence of the exqtansion of the mmber density function nj^Cz.t) in 

-erms of orthogonal polynomials with a given weighting functions in the range (I,s) Css2) is 

sufficiently npid and for droplet mass grid points ^V4.i*4x. (s$2) the approximation of nj^(z,t) 

by means of the first two terms of the expansion (Ij *is ^sufficiently correct. 
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